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ANALYTICAL STUDY OF DIFFUSION -CONTROLLED CHAR OXIDATION 
AND ITS EFFECT ON STEADY-STATE ABLATION 


OF PLASTIC MATERIALS 

By Gerald D. Walberg 
Langley Research Center 

SUMMARY 'X \^>\ 

Z) 

An analysis of diffusion-controlled char oxidation has been carried out and differs 
from previous analyses in one important respect: previous investigators have assumed 
that only the carbon present in the char would react, but the present analysis has shown 
that the carbon present in both the char and the pyrolysis gases (i.e., all the carbon orig- 
inally present in the virgin plastic) reacts and that the reaction rate is determined by 
the rate at which oxygen can be supplied by diffusion through the boundary layer plus 
pyrolysis of the plastic. Because of this fact, the char oxidation rate, the steady-state 
ablation rate, the boundary -layer composition at the char surface, and the heat released 
by chemical reactions were found to be dependent on the composition of the virgin plastic 
and the pyrolysis gases. 

For a given flow environment, the char oxidation rate and the steady -state pyrolysis 
rate were found to increase as the hydrogen, nitrogen, and oxygen contents of the virgin 
plastic increased and to decrease as the carbon content of the virgin plastic increased. 

It was found that the heat released by both gas phase and heterogeneous reactions can be 
accounted for by an "effective enthalpy potential" which is a function of the stream oxygen 
mass fraction and the compositions of the virgin plastic and the pyrolysis gases. Partic- 
ularly high heat releases are predicted for polymers which contain nearly equal numbers 
of carbon and oxygen atoms and produce pyrolysis gases containing large quantities of 
C2H2 and C2H4. 


INTRODUCTION 

Since the efficiency of a charring ablator is determined to a large extent by the char 
removal rate, an understanding of char removal mechanisms is essential for the predic- 
tion of the performance of a charring ablator. For the steady -state ablation case, the 
linear rate of pyrolysis is equal to the char removal rate; hence, the virgin plastic 
pyrolysis rate is known if the char removal rate can be determined. 


When the char layer under consideration possesses adequate structural strength so 
that mechanical failure (removal by aerodynamic shear, spallation, etc.) does not occur, 
the char removal rate is determined by the rate of chemical reaction between the char, 
the gaseous pyrolysis products, and the chemical species present in the boundary layer. 
Since an ablator which loses large amounts of char due to mechanical failure is not suit- 
able for use on reentry heat shields, it is not unreasonable to expect that chemical 
reactions will be the primary char removal mechanism during most reentries. 

The maximum rate at which char oxidation can proceed is ultimately determined 
by the rate at which fresh supplies of oxygen can be transported to the char surface. 

When a laminar boundary layer exists over the char surface, the transport mechanism of 
interest is molecular diffusion and this maximum rate is referred to as the diffusion- 
controlled oxidation rate. Since carbonaceous materials possess extremely high intrin- 
sic chemical reactivities at temperatures above 3000° R (1667° K) (refs. 1 and 2), the 
char-layer oxidation is generally assumed to be diffusion controlled under most reentry 
conditions. (See refs. 1 to 5.) 

In the present report, an analytical study of the diffusion-controlled oxidation of an 
ablative char layer is presented and, within a range of pressure and temperature which 
is typical of stagnation-point conditions during orbital reentry, an approximate closed- 
form solution is obtained. Although this solution is based on a number of simplifying 
assumptions, it is believed to be a sound conceptual picture of the diffusion-controlled 
oxidation process which demonstrates the role of pyrolysis gases in determining the 
oxidation rate (and hence the steady-state pyrolysis rate) and provides a relation between 
the original chemical composition and the cold-wall effective heat of ablation of a 
charring ablator. 


SYMBOLS 

All dimensional quantities are presented in terms of both U.S. Customary Units and 
the International (SI) System of Units. (See ref. 6.) 

A,B,C,D stoichiometric coefficients per mole of boundary -layer gas mixture 

(see eq. (100)) 

a,b, . . .,o equilibrium coefficients (see eqs. (71) to (85)) 

Cp specific heat of boundary -layer gas mixture, British thermal units per 

pound-°R (joules per kilogram-°K) 
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D 


coefficient of diffusion, feet^ per second (meter per second) 
mass of char produced by decomposition of unit mass of virgin plastic 


f 

H 

AH 

AH r 

AH f>i 

h 

h m 

Ki 

Ki 

Ki 

K P,(n) 

k 

l 

M 

a 

M 


stream enthalpy, British thermal units per pound (joules per kilogram) 

enthalpy potential across boundary layer, H e - H w , British thermal units 
per pound (joules per kilogram) 

overall heat of reaction, British thermal units per mole of carbon (solid) 
(joules per mole of carbon (solid)) (see eqs. (113) and (116)) 

heat of formation of species i, British thermal units per mole of i (joules 
per mole of i) 

convective heat-transfer coefficient, pounds per foot^- second (kilograms 
per meter^-second) 

mass-transfer coefficient, pounds per foot^ -second (kilograms per 
meter^- second) 

mass fraction of species i, mass of i per unit mass of gas mixture 

elemental mass fraction, mass of chemical element i (no matter what 
chemical form it may appear in) per unit mass of gas mixture 

mass of species i per unit mass of gases injected into boundary layer 

equilibrium constant (expressed in terms of partial pressures) for 
reaction presented in equation (n) 

thermal conductivity of boundary -layer gas mixture, British thermal units 
per foot- second -°F (joules per meter-second-°K) 

exponent in transpiration cooling expression (see eqs. (125) and (127)) 

molecular weight of gas mixture at char surface 

average molecular weight of gases injected into boundary layer 


3 



V 


Mi 


m 


m C 

m H 

m N 


m Q 


N 


A 


N 


Le 


n ik 


P 


Pi 


Q* 


Qeff 


q 


q cw 


^net 


molecular weight of species i 

mass rate of oxidation or ablation, pounds per foot^ -second (kilograms 
per meter2 -second) 

mass of carbon per unit mass of virgin plastic 
mass of hydrogen per unit mass of virgin plastic 
mass of nitrogen per unit mass of virgin plastic 
mass of oxygen per unit mass of virgin plastic 


Avogadro's number 

PDf2 c p 

Lewis number, — : 

k 

number of moles of species i per mole of species k (see eq. (23)) 

pressure, pounds per foot^ or atmospheres (newtons per meter2) 

partial pressure of species i, pounds per foot^ (newtons per meter^) 

thermochemical heat of ablation, British thermal units per pound (joules 
per kilogram) (see eq. (142)) 

q 0 

cold -wall effective heat of ablation, - , cw , British thermal units per 
pound (joules per kilogram) 

convective aerodynamic heat flux, British thermal units per foot^ -second 
(joules per mete r^- second) 

cold-wall aerodynamic heat flux, British thermal units per foot^-second 
(joules per meter^ -second) 

net heat flux to surface of ablator, used in defining Q*, British thermal 
units per foot^-second (joules per meter^ -second) (see eq. (142)) 
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heat flux due to both convection and chemical reactions, British thermal 
units per foot^- second (joules per mete r^- second) 

r nose radius of hemispherical body, feet (meters) 

T temperature, °R (°K) 

v boundary -layer velocity normal to char surface, feet per second (meters 

per second) 

W ratio of diffusive oxygen flux to char oxidation rate, defined by 

equation (95) 

Xi mole fraction of species i, moles of i per mole of gas mixture 

Xj moles of species i per mole of pyrolysis gas 

y distance normal to char surface, positive when directed away from char, 

feet (meters) 

rj transpiration cooling factor (see eq. (120)) 

A ratio of virgin plastic pyrolysis rate to char oxidation rate (see eq. (7)) 

p boundary -layer density, pounds per foot^ (kilograms per meter 3) 

a,T,v,(p indices in empirical polymer formulas, indicating the respective number 

of carbon, hydrogen, nitrogen, and oxygen atoms in the plastic 


Subscripts: 

C,H,N,0 pertaining to carbon, hydrogen, nitrogen, and oxygen, respectively 

c pertaining to ablative char layer 

ch pertaining to or resulting from chemical reactions 

e evaluated at outer edge of boundary layer 
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i,j,k pertaining to ith, jth, or kth chemical species or element 

s evaluated at stagnation conditions 

vp pertaining to virgin plastic 

w evaluated at char surface 

12 denoting values for binary gas mixture 


Superscripts: 

o evaluated for zero mass injection into boundary layer 

* steady -state solution 


ANALYSIS 


In this section, the governing equa- 
tions and the assumptions used in the pres- 
ent analysis are stated, and the ranges of 
temperature and pressure within which an 
approximate closed-form solution is 
obtainable are identified. The distance y 
used in this analysis and the various zones 


Mass fraction 
boundary j 
layer I 

Char 


Virgin 
plastic [ 


Boundary-layer Reaction 
y species ^ /-'"products 



-Outer edge of 
boundary layer, e 


^/~Char surface (w) 
/-Pyrolysis interface 


L-64- 10,940.1 

Figure 1.- Charring ablator showing zones and interfaces. 

and interfaces pertinent to the behavior of a charring ablator as analyzed herein are pre- 
sented in figure 1. The virgin plastic is assumed to contain only nitrogen, oxygen, carbon, 
and hydrogen. 


Governing Equations 

Consider the decomposition of a hydrocarbon polymer according to the equation 

- Pyrolysis gas + C(s) (1) 

At the char surface, mass conservation of the element nitrogen may be represented by 
the equation 
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-h 


hi,n( : 


K 


N,e " ^N,w 



Diffusion away from 
char surface 


Convection 
away from 
char surface 


vM 


N 


+ tM h + + <pMqJ 


m. 


vp 


( 2 ) 


Flow of pyrolysis gases 
from pyrolysis zone to 
char surface 


where o, t , v, and cp represent the relative number of atoms of carbon, hydrogen, 
nitrogen, and oxygen in the virgin plastic, and h m is the elemental mass-transfer 
coefficient for nitrogen, which is defined by the equation 

( pD N ~of) w = h m,N( K N,e " K N,w) (3) 

The quantity in parentheses on the right-hand side of equation (2) represents the mass of 
elemental nitrogen contained in the pyrolysis gases produced by decomposition of a unit 
mass of virgin plastic. 

In order for equation (3) to be rigorously valid, all the species containing nitrogen 
must have nearly equal diffusion coefficients. Since subsequent equations will place 
similar requirements on the species containing oxygen, hydrogen, and carbon, the use of 
the diffusion coefficients Djj, Dq, Djj, and D<-. implies that the diffusion coefficients 
of all the species are nearly equal - that is, the gas mixture is effectively binary. This 
binary mixture requirement is discussed in the section "Review of Assumptions." Equa- 
tions similar to equation (2) may be written for oxygen, carbon, and hydrogen as follows: 


-Vo(Ko,e - Ko,w) - H- ' 


(4) 


-h m ,c(Kc,e - Kc,w) + (*CP V ) W = 


a M, 


loM c + tM h + + <pM Q 


- f irhyp + (5) 


-V,H(K H ,e - K„, w ) + M, - 


(pM. 


m 


vp 


( 6 ) 


where f represents the mass of char produced by decomposition of a unit mass of 
virgin plastic. Now define A as 


m vp 

m c 


(7) 
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and note that 


(pv) w = m c [l + A(1 - f)] 


( 8 ) 


The respective masses of nitrogen, oxygen, carbon, and hydrogen per unit mass of virgin 
plastic can be expressed as 


I'M 


m N = 


N 


oM^, + tMjj + I'M^ + c^Mq 


(9) 


<P M, 


m Q = 


O 


aM c + tM r + i4V1 n + c/9Mq 


( 10 ) 


ctM, 


mn = 


+ tM h + + cpM Q 


(ID 


tM. 


m H = 


H 


ctMq + tMj^ + I'M^ + cpM Q 


(12) 


Substitution of equations (7) to (12) into equations (2), (4), (5), and (6) yields the following 
expressions for the elemental mass fractions at the char surface ^by definition, 

^C,e = %,e 


= 0 ): 


S- k N,e h m,N + m N Xtil c 

Kp. h„ ~ + m n Ara„ 

Jr _ Q,e m,Q U c 

’ W h m,0 + T 1 + X(1 " f) ] “c 



n 

1 

u 

.a. 


m 

c 

h p + 

m j C 

[l + X(1 

- 

3 

ril c 


_ m^Am c 

H ’ W h m,H + C 1 + X ^ " f )] m c 


(13) 


(14) 


(15) 


(16) 
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Equations (13) to (16) are the governing equations for this study, 
however, cannot be solved unless the mass-transfer coefficients h m> j 


These equations, 
can be evaluated. 


Evaluation of Mass-Transfer Coefficient 

In this report, the mass -transfer coefficients h m j are expressed in terms of 
the convective heat -transfer coefficient through the application of an analogy between the 
transfer of energy and mass in a laminar boundary layer. This analogy (with minor 
variations) is often used in treating the diffusion-controlled oxidation of graphite. (See 
refs. 2, 3, 7, 8, and 9.) 

Lees in reference 9 shows that, for the case of a laminar, unit Lewis number, 
stagnation -point boundary layer, the ith elemental mass fraction and the total enthalpy 
are related by the following equation: 


Hence 


and since 


^i “ ^i,w _ H s - H s w 


Ki,e - Ki w H s e - H S;W 


Kj.e - *i.w 
9y "H s e -H s w 3y 


PD^ =1 
Le k 


(17) 


(18) 


(19) 


it follows that 


9Kj /~ ~ \ 

pD 12 ( K i,e ' K i,w) N 


q 


9 y 


Le 


Hs,e " Hs,w 


" (^i,e ^i,w)k 


But the ith elemental mass-transfer coefficient is defined by the equation 


( 20 ) 


9K; /- ~ X 

pD i - ( K i,e ' KijwJbm,! 


9 y 


( 21 ) 


Hence, if the boundary -layer mixture is effectively binary, 


^m,i “ k 


( 22 ) 


This assumption is discussed in the section "Review of Assumptions.’ 
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Pyrolysis Gas Composition, Char Composition, Assumed Reactions, 
and Equilibrium Expressions 

The elemental mass fractions Kj which appear in equations (13) to (16) can be 
expressed in terms of the mole fractions of the various chemical species by the following 
equation: 


Ki = 



where 


(23) 


Mj atomic weight of elemental species i 
M molecular weight of gas mixture 

Xk mole fraction of kth chemical species in gas mixture 
nii, number of moles of i per mole of k 


The chemical species which must be considered are determined by the reactions which 
take place between the pyrolysis gases, the boundary -layer species, and the char. 

The elemental mass fractions K^ w are expressed in terms of the reactants and 
products under the assumption that chemical equilibrium (both heterogeneous and gas 
phase) exists at the char surface. The pyrolysis gases are assumed to be composed of 
Ng, Hg, NHg, CO, COg, HgO, C 2 H 2 , C 2 H 4 , and CH 4 , and the char is assumed to be 
graphite. 

The following chemical reactions are considered: 


2NH 3 z 3H 2 + N 2 

(24) 

CH 4 + 20 2 z C0 2 + 2H 2 0 

(25) 

C 2 H 2 + |°2 = 2 C0 2 + H 2° 

(26) 

C 2 H 4 + 20 2 Z 2CO + 2H 2 0 

(27) 
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(28) 


f°2=° 


£ H 2= H 

(29) 

|h 2 + |o 2 zoh 

(30) 

H 2 + |°2= H 2° 

(31) 

|n 2 =n 

(32) 

iN 2+ |o 2 =NO 

(33) 

c(s) r c(g) 

(34) 

c(s) + 1 o 2 r co 

(35) 

C(s) + C0 2 z 2CO 

(36) 

C(s) + h 2 o -h 2 + co 

(37) 

c(s) + n:cn 

(38) 

C(s) + NH 3 r HCN + H 2 

(39) 


The notation s or g indicates species in a solid or gaseous state, respectively. The 
equilibrium constants for reactions (eqs. (24) to (39)) may be written as follows: 


*S>,(24) - 


% ^2 X] 

PNHo 2 




H 2 N 2 2 
x nh 3 2 


(40) 


K 


PcOoPHoO X C0„ X H„0 


'2 2 


2 2 


p ’ (25r Pch 4 Po 2 2 


x ch 4 V 


(41) 
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(42) 


P C0 2 2p H 2 Q _ X CQ 2 2x H 2 Q _ 1/2 

P>(26) " p C H p O 5/2 X C H X 0 5/2 P 
^2 rt 2 u 2 °2 n 2 u 2 


K. 


P? (27) 


p C0 2p H 2 0 2 X C0 2 x H 2 0 2 

Pc 2 H 4 P °2 2 Xc 2 H 4 X °2 2 


(43) 


Pr 


0 1/2 


Kp ,(28) = -^ = 7-i^ p 


p o 2 - -o 2 


(44) 


K- - P H _ X H _ 1/2 

P> (29) - p 1/2 - x 1/2 P 

H 2 H 2 


(45) 


K, 


p OH 


OH 


P ^ 30 ^ P H 1//2 P n X„ 1 / 2 X n x / 2 

h 2 u 2 h 2 u 2 


(46) 


P H 2 ° 


X T 


K - « - H 2° p -l/2 

K p,(31) i“7o i /o p 


P H 2 P 0 2 1/2 ~ X H 2 V /2 


(47) 


K _ P N _ X N 1/2 

p,(32> "pn 2 1/2 " x n 2 1/2 


(48) 


K 


PNO 


X 


NO 


P,(33) ' p N ^l/2p 0 ^l/2 " 5^172^172 


(49) 
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K p, (34) ~ P C(g) “ X C(g) p 


(50) 


if _ P CO _ X CO 1/2 

Kp ’^ 35) Po ^ 2 x o 1//2 
u 2 U 2 


r = P CQ 2 _ X CQ 2 

p ,( 36 ) " p C o 9 _ x COo 


PHgPCQ X H 2 X CO 


“^>-*55-^ 


P CN X CN 


P,(38) Pj^ x N 


^P»(39) ^NH 


p hcn p h 2 x hcn x h 2 


Equations (40) to (55) may now be used to express the mole fraction of each species in 
terms of the mole fractions of N 2 , H 2 , and CO. The resulting expressions are 


X — aX 3 / 3 y 
X NH 3 - aX H 2 X N 2 


X CH 4 " bX H 2 


Xp TT = eXtt 
u 2 M 2 H 2 
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x r H = dX„ ' 
C 2 H 4 H 2 


X Q " eX CO 


X — fX 
X H ~ 1X H 2 


X OH “ gX CO X H, 


1/2 


X H 2 0 " hX H 2 X CO 


X - iX ^/^ 
A N " lx n 2 


X NO " ^ X N 2 1/2x CO 


X c(g) " k 


x o 2 =* x co 


x co 2 " mX CO 


X CN " nX N 2 ^ 2 


X HCN - oX H 2 1/2x N 2 V2 


(59) 

(60) 
(61) 

(62) 

(63) 

(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 


(70) 


where 


a = 


P_ 

K P,(24) 


1/2 


K P>(35) 4 P 

K p,(37) 2K p,(36) K p,(25) 

K 5 

P> (35) 

K P,(36) 2k p,(37) K P ,(26) 
K- 4 

. K P,(35) P 

d 2 

K P» (27) K p, (37) 

c _ K P,(28) 

*^,( 35 ) 


1 * K P,(29)P' V2 
K P,(35) 


(71) 


(72) 


(73) 


(74) 


(75) 


(76) 


(77) 


h = 


K P,(37) 


i=K p,(32)P 


- 1/2 


(78) 


(79) 
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( 80 ) 


*»,<»» ) pV2 

K p,(35) 


k = K. 


P,(34)P 




(36) 


n-K p,(38) K p,(32)P' 


P>(39) 
S, (24, 1/2 


The governing equations (eqs. (13) to (16)) may now be restated as follows: 


= 3/2 X m X / 2 + jX XT 1/2 X™ + nX, T 2 + iX XT V 2 


H 2 N 2 


j X N 2 X CO + nX N 2 


+ oX h 2 V2x n 2 1/2 


K N,e h + m N Xlh c 
h + [l + A(1 - f)Jm ( 


£ °' w = M ( JXn 2 1/2Xc0 + 2iXc ° 2 + 2mXc ° 2 + eXc ° + ^CX^h/ 72 


. kv v . v \ K 0,e h + m O A “c 

+ hx H« x co + x co ~ r f 77: — 7 rnr~ 
2 J h + [1 + X (1 - f)J m 
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K C ,w = ^( bX H 2 2 + Xco + 2CX H 2 * k + "*C0 2 + ^2 V2 + <*H 2 1/2 *N 2 V2 


+ 2dX 


H„ 


'w 


I + 

( m c - f ) 


“c 

h+ [ 

1 + A(1 - 

f) 

]™C 


( 88 ) 


K„ ... = ^3aX u 3/2 X„ V2 + 4bX H2 2 + 2X H2 ♦ 203^ + fX^ 1 / 2 + gX^^ 2 


H.«' = mT H 2 N 2 


+ 2hX„ X™ + oXtt 1//2 X xt 1//2 + 4dX 2 


"h 2 ~co t u ' h 2 ~n 2 


Hr 




m H Xm c 


A(1 - f j] m c 


(89) 


Approximate Solution 

General expressions .- When all the equilibrium coefficients (a, b, etc.) are «1.0, 
equations (86) to (89) may be approximated as follows: 


- ,^n x l%l,e h * m N Xm c 

N,» S N 2 ,w N 2 > w ~h + + A(1 -f)]m. 


(90) 


M M„ 

V"J X CO,w K 


KQ,e h + V A C 

M CO € °’ W h + t 1 + M 1 - f)] m f 


(91) 


M„ 

Kp ... = X, 


M, 


K, 


C 1+ 

( m c - f ) A ] 

“c 

h + | 

1 + A(1 - f) 

]“c 


(92) 


2 M 


K 


H 


X„ ... = K, 


m„Am„ 
H c 


H,w - jj H 2 ,w H 2> w - h + Q + X(1 _ f) ] A 


(93) 
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and equations (91) and (92) may be solved 

“c 

where 

W = — f + ( m c - f ) A ] - m o A {95) 

c 

It should be noted that the solution for m c , presented in equations (94) and (95), is valid 
for the general case of transient ablation. The parameter W is, however, a function 
of X and in the transient case X is not known but must be obtained from experiment 
or through the iterative solution of a transient ablation computer program. For the 
steady-state ablation case, however, X is a known constant. The steady-state ablation 
case is considered in detail in the subsequent section of this paper. 

The boundary -layer composition at the char surface can be obtained from equa- 
tions (90), (91), and (93) as follows: 


for m as follows: 

V 


hK o,< 


w 


( 94 ) 


K %w 


K CO,w “ 


K H 2 ,w = 


W^ + m,X 

Ko.e N 


w + 

1 + X(1 - 

f) 

Ve 



M CO 

jw + itIq) 


M 0 

3 

W + 

1 + X(1 - 

f) 

Ve 




m H X 



W 


K 


+ 1 + X(1 - f) 


0,e 


(96) 


(97) 


(98) 


The mass fractions IC can, of course, be expressed in terms of mole frac- 

i,w 

tions X, through the relation 
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( 99 ) 


M; 


X* - L- 

Al - K; 


y 3. 

Z,Mj 


] 

Once the boundary -layer composition at the char surface has been determined, an 
overall heat of reaction may be calculated as follows. An equation for the overall reac- 
tion per mole of boundary-layer gas mixture may be written as 


[ao 2 + BN 2 ] + 


X CO, X CH. X C,H, X C 0 H, 

c (CO +^— £ co 2 + - CH 4 + _ Z z C 2 H 2 + 4 c 2 h 4 

x co x co x co x co 


X N, X NH„ X H„0 X H 0 

+ =-^ N 9 + — — £ NH, + H„0 + = -^- H 0 

Y ^ y y m y m 

A co A co A co A oo 


[DC(s)] ; [x N2jW n 2 + x COjW co + x H2iW h 2 J 


(100) 


where is the number of moles of species i per mole of pyrolysis gas. The terms 
of equation (100) are identified as follows: 


^Species from stream] + [Pyrolysis gases] + [char] ~ [Reaction products] 


The requirement that the number of atoms of a given element must be the same on both 
sides of equation (100) results in the following relations: 


2A + 


C 1 + 


X CO x co; 


= X 


CO,w 


( 101 ) 


2B 



= 2X. 


N 2 ,w 


( 102 ) 
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/ x co 2 x ch 4 2x c 2 h 2 2x c 2 h 4 \ 

C 1 + — — - + i + _ 1 - + ■ _ 2 4 + D = X 

\ x co x co x co x co / 


CO,w 


*CH 

2x r h 
C 2 H 2 

4x r h 
C 2 H 4 , 

3X NH 3 

— + 


+ — — + 


\ x co 

x co 

X 

o 

o 

x co 


2x h,o 2x h, 
+ _ ■■■ " + — — - 

x co x co, 


H 2 .w 


but, for each mole of reaction products produced, N^C 11 + — + 


2x co 2 x h 2 o 


x co x co 


oxygen 


atoms entered the reaction in the form of pyrolysis gases and N^D atoms of carbor 
entered the reaction in the form of char. Hence, from the definition of f and mQ, 


\ x co x CO / _ m O M C x 

r» Mm 


Similarly, 


C 1 + ^ “ + • • 

\ 


= ( m c - f ) x 


c — — - + 
\ Xmm 


m H M C x 


/ 2x n 

cl — — + 

i CO 


m N^C x 
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and equations (101), (102), (103), and (104) with the aid of equations (105) to (108) can be 
expressed per mole of C(s) (char) by dividing through by D. Hence, 


2 A . m O M C x .. X CO,w 

D Mq D 


(109) 


2X 

2B . m N M C X _ N 2 >w 
D ' M n d 


( 110 ) 



, 1 - X CO.w 


( 111 ) 


2x 

hirMcX ®2 ,w 

M H D 


( 112 ) 


Since the overall heat of reaction is defined as 

AH^ = Enthalpy of products - Enthalpy of reactants 

from equation (100), the overall heat of reaction per mole of solid carbon (char) consumed 
may be stated as 


AH R " C q’ W AH f,CO ‘ DX^O^CO + X C0 2 AH f,C0 2 + X CH 4 AH f,CH 4 
+ X C 2 H 2 AH f,C 2 H 2 + X C 2 H 4 AH f,C 2 H 4 + X NH 3 A ®f,NH 3 + X H z O AH f,H 2 oj 


(113) 


but according to equations (111) and (105) 
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‘ C lP = (m C - *)* + 1 


(114) 


c m O M C x 


DX C0 M 0 X CQ + 2X C02 + X H20 


(115) 


Hence equation (113) becomes 


AH r = |m c - f)x + lj AHf co 


m O M C x 


2 s l AH f,i 


° x co + 2x co 2 + x h 2 o 


The heat flux due to chemical reactions (recombination reactions for oxygen and nitrogen 
are not included since they are usually accounted for in computing convective heat fluxes) 
may now be expressed as 


^ch _ 


AH R m c 


hKo,e 

M C W 


Y, x i AH f i 

A, m O M C x i 1 M 

r c ■ f)x + r 1 ^ 0 • -®5- w 2 x„, ,x 


(117) 


co T co 2 T h 2 o 


(The minus sign is used in equation (117) since, by definition, an exothermic reaction has 
a negative heat of reaction.) Equation (117) may now be used to define an "effective 
chemical enthalpy potential AH c ^" as follows: 


AH ch 


~ f 1 

_ q ch_ K 0,ek f a m O M C x i 

- — - - M^wf “C - ')* * IJAH^CO - J—W— 


co + 2x co 2 + x h 2 o 


(118) 


Hence, the total heat flux (convective + chemical) is expressible as 


q t = h(AH + AH ch ) 


(119) 
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Since expressions for the gas composition at the char surface have been derived, 
consideration of the transpiration cooling effect produced by the pyrolysis gases is of 
interest. In practical applications the reduction in heat transfer which results from the 
injection of gas into the boundary layer is often accounted for by using the relation 


h = h° - ?7(pv) w 


( 120 ) 


and in reference 10 it is suggested that the transpiration factor rj varies approximately 

* 1 A 

as M where M is the molecular weight of the injected gas and l ranges from 
1/3 to 1/4. If the gas phase combustion reactions all occur near the char surface, the 
species which are actually injected into the boundary layer are the combustion products, 
that is, Ng, Hg, and CO. The mass fraction of N 2 in the injected gases is expressible as 


and similarly 



m 


N 0 


m XT ni 
N vp 


m N X 


(pv) w (1 - f)m yp + m c 


(1 - f)A + 1 



m »2 _ m H™vp _ m H X 
(pv) w (1 - f)m yp + m c (1 - f)A + 1 


( 121 ) 


( 122 ) 
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CO 


_ m co _ ( 


m O + m c - f)A vp + rn ( 


(pv) 


w 


(* - f )* V p +,h c 


_( 


m O+m c 




A + 1 


(1 - f)A + 1 


(123) 


Hence, the average molecular weight of the injected gases is given by the expression 


M = ( 1 l i ) - x l - 1 (124) 

m^jA m jjA (iHq + - f^A + 1 

M N 2 m h 2 m co 


Now the ratio of the transpiration cooling factor for a gas of arbitrary molecular weight 
to the transpiration cooling factor for carbon monoxide (M = 28) is 


23 


V 


t? M= 28 


(mV' = 

(1 - f)X + 1 

-l 

13mTjX 
1 i H 

\28/ 

m^X + 14mjjX + m^X + m^X - fX + 1 


(1 - f)X + 1 


( 125 ) 


and since for a laminar boundary layer 


?7 M=28 ~ ^air 


0.7 


(126) 


equation (120) may be approximated by the following expression 

0.7 |1 + S—\ (p V ) 




h ~ h 


(1 - f)X + 1 


'w 


(127) 


Steady-state expressions .- The solutions obtained in this analysis (as presented in 
eqs. (94) to (98), (116) to (119), and (127)) are applicable to the case of transient ablation - 
that is, the case where the char surface and the pyrolysis interface move at different 
linear rates. However, these solutions depend on the ratio of pyrolysis rate to char 
oxidation rate X, and X is not known a priori except in the case of steady-state abla- 
tion where X = 1/f. The limiting case of steady ablation is, however, an important one 
and hence the solutions for X = 1/f are presented in the following equations: 



m vp = 


*0.e h 


fW* 


(128) 


where 


* M o 

fW* = m^, - m~ = 


M n (o - <p) 


M, 


C “O 


oMq + tM h + + <pM Q 


(129) 




( 130 ) 
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Kh 2-” ’ ^ ^ Kc °.» = ™ Kc °.»' 


(131) 


K, 


>‘+"o) 


CO, w fw* 

~-+ 1 

*0,6 


(132) 


m. 


m. 


AH n = — ^ AH, 


O 3 


1 


R f f,C0 f 4\x 2X + X 

1A co + za co 2 + A H 2 0/ 


(133) 


AH, =- tt = 

ch h 


| m 0 


I * 1^4 


(f m c- m o) 


m 0 AH. 


■c f rn — — — — 

191’ m m i\ ’ X rn + 2X r , r , +X„ n 
12| — m c -m 0]\ CO C0 2 HgOj 


J 5 a< 


/ 


16(cr - cp) 


x i aH f,i 


CO ” — — — 

V ’ x co + 2x co 2 + x H 2 o, 


(134) 


h = h° - rjm vp = h° - 0.7 ^1 + 13m H ^ m yp (135) 

DISCUSSION 

In this section, the assumptions used in the analysis are reviewed, a physical 
interpretation of the approximate solution is presented, and the effects of virgin plastic 
composition and pyrolysis gas composition on the char oxidation rate, the steady-state 
ablation rate, the cold-wall effective heat of ablation (defined subsequently), and the heat 
release due to chemical reactions are discussed. In the section "Analysis," both general 
and steady-state expressions of the approximate solution were presented. Before a 
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solution to the general expressions can be obtained, the ratio of pyrolysis rate to char 
oxidation rate X (which is not known except for the steady-state ablation case) must be 
specified. Hence, while the general expressions can be used in conjunction with a com- 
puter program (such as that presented in ref. 11) to yield transient ablation rates, they 
do not lend themselves to demonstrating the general characteristics of diffusion- 
controlled char oxidation. For steady -state ablation, on the other hand, X = 1/f and as 
a result the equations are simple and provide a clear physical picture of diffusion- 
controlled char oxidation. Accordingly, the following discussion is restricted to the 
steady-state case. Any remarks concerning the char oxidation rate apply equally well 

to the virgin -plastic ablation rate since m = fm . 

c v p 

Review of Assumptions 

The major assumptions used in the present analysis are: 

(a) The virgin material is composed of only oxygen, hydrogen, nitrogen, and 
carbon. 

(b) The mass -transfer coefficients for all gas phase species are numerically equal 
to the convective heat -transfer coefficient. 

(c) All the equilibrium constants in equations (86) to (89) are «1.0. 

(d) Chemical equilibrium (both gas phase and heterogeneous) is obtained at the char 
surface . 

(e) The transpiration factor r\ varies inversely with the kth power of the mean 
molecular weight of the gases injected into the boundary layer. 

(f) The pyrolysis gases are composed of Ng, Hg, NHg, CO, COg, CgH^, CgHg, 

CH^, and HgO. 

Most presently used high -temperature resins (silicone resins and rubber are 
notable exceptions) are composed of O, H, and C. Nylon, which is commonly used as a 
reinforcement in ablative materials, is composed of O, H, N, and C; hence, mixtures 
such as epoxy nylon or phenolic nylon can be treated by the techniques presented herein. 
Also, the present method could be applied to a material composed of a hydrocarbon resin 
and an inert filler. 

The assumption made in evaluating the mass -transfer coefficients is commonly 
made in analyses of carbon combustion. (See refs. 2, 3, 7, 8, and 9.) In reference 9 
this assumption is shown to be valid whenever the gas mixture can be considered binary 
with respect to diffusion and the Lewis number is unity. Even though neither of these 
requirements are rigorously met for carbon combustion, use of the mass -transfer - 
coefficient assumption yields results which are in qualitative agreement with experiment. 
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A measure of the validity of this assumption for carbon combustion can be obtained 
by comparing the present solution with that presented in references 1 and 4. The 
analysis presented in references 1 and 4 considered multicomponent diffusion and placed 
no restriction on Lewis number. From that analysis the following expression for the 
diffusion-controlled stagnation -point carbon oxidation rate in air was obtained and pre- 
sented in reference 4: 

= 6.35 x 10' 3 lb/ft 3 / 2 -sec -atm 1 / 2 = 0.05378 x 10' 3 kg/m 1//2 -N 1//2 -sec (136) 

Before the results of the present analysis can be compared with equation (136) a value 
for the heat -transfer coefficient must be obtained. Evaluation of the Fay and Riddell 
relation (ref. 12) has shown that values on the order of 

h°^^= 3.92 x 10" 2 lb/ft 3 / 2 -sec -atm 1 / 2 = 0.033198 X 10 -2 kg/m 1 / 2 -N 1 / 2 -sec (137) 

are typical of stagnation-point conditions on a manned vehicle during entry from earth 
orbit. If this value and the oxygen mass fraction for air (0.232) are substituted into 
equations (94), (95), and (127) there results, 

= 6.07 x 10" 3 lb/ft 3 / 2 -sec -atm 1 / 2 = 0.051407 x 10" 3 kg/m 1 / 2 -N 1 / 2 - sec (138) 

which agrees with equation (136) to within 5 percent. The total heat flux calculated from 
equation (119) was also found to agree with the values presented in references 1 and 4 
to within 5 percent. Hence the unity Lewis number assumption may be considered valid 
for carbon combustion. Only the presence of light gases (primarily H 2 ) makes the 
assumption less applicable to the present analysis than to carbon combustion. Hence the 
assumption should be quite good for low values of free -stream oxygen mass fraction 
^as *0,e- 0 and Kg w — 0, eq. (98) j and although the validity of the assumption 

decreases with increasing *0,e> it should provide an acceptable qualitative description 
of the transport process for air. This assumption is considered in the light of calculated 
char-surface boundary -layer compositions subsequently. 

The approximate solution was obtained by assuming that the equilibrium coefficients 
appearing in equations (86) to (89) were small compared with unity. A more rigorous 
statement of this requirement would be that the sum of all the equilibrium coefficients 
appearing in a given equation must be small compared with unity. In order to assess the 
validity of this assumption, the equilibrium coefficients were evaluated by using the data 
presented in reference 16. The variation of these coefficients (evaluated for a pressure 
of 1 atm) with temperature is presented in figure 2. By using the curves presented in 
figure 2, ranges of char surface temperature and pressure for which the sum of the equi- 
librium coefficients appearing in any one of equations (86) to (89) was <0.1 were identified. 
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The ranges of char surface tem- 
perature and pressure within which 
the sum of these coefficients is 


d (c 2 h 4 ) 

rb(CH. ) 


<0.1 (and, hence, within which the 

i n (CN) 

approximate solution is valid) are ' — / / /“ f(H) 

presented in figure 3. The char — * , l r c(9 u )) 

/r C ( c 2 ) 

surface temperatures and pres- o 1 

sures experienced by a manned syl' / 

ballistic vehicle (such as the I / / /f / 

Mercury spacecraft) during entry ^ / / // / 

from earth orbit fall within this \ / /// / /' lN 

range and hence the present " 2— \V / // / / / 

approximate solution should be 1 r ySJ jj / / 

applicable to such vehicles. On § _ 3 — Vv\ < / 

the other hand, ballistic missiles | / 

and lunar return vehicles experi- 1 4 V yw\ / / 

ence temperatures higher than /\\j\ d(c h ) 

those shown in the cross-hatched ^ 

region of figure 3. Solution of the ~ 5 / /^Xg ( ( °H) 

governing equations at these r - I / j [f\ 

higher temperatures is straight- -e— 1/ ! 

forward but would require the use — / i I / / J (N0) 

of a digital computer. _ 7 / // ^ m(co 2 ) 

The assumption of chemical — / 

equilibrium at the char surface is ° ■ i ' l 1 I I J I i I ! I 

, „ 0 1000 2000 3000 4000 5000 6000 °K 

a natural consequence of con- ( ^ ^ ( ^ 

sidering the case of diffusion- 0 1800 3600 5400 7200 9000 ° R 

Temperature , T 

Controlled Char oxidation. In any Figure 2.- Variation with temperature of equilibrium coefficients, 

kinetic system, the slowest step Pressure, 1 atm (0.1013 MN/m2). 

is rate controlling. Hence, in the diffusion-controlled case, the rate of transport of 
reactants must be slow compared to the rates of both homogeneous and heterogeneous 
reactions. Since the time required for the slowest reaction to equilibrate is small com- 
pared with the time required to deliver fresh quantities of reactants to the char surface, 
chemical equilibrium is obtained. 


Temperature , T 

Figure 2.- Variation with temperature of equilibrium coefficients. 
Pressure, 1 atm (0.1013 MN/m2). 


The assumption made with regard to the transpiration factor tj is semiempirical 
in nature and is based on experiments and analyses involving nonreacting boundary 
layers. This assumption would be applicable to the present problem if all the gas phase 


28 


reactions occurred in a N 

thin zone near the char -, wZ *1™ 

surface. When this sit- 
uation occurs, the gases 
injected into the bound- 
ary layer are the reac- 
tion products which 
diffuse from the reac- 
tion zone through the 
boundary layer without 
undergoing further reac- 
tion. It should be noted 
that, in analyzing reentry 
heating, several investi- 
gators have assumed a 
frozen boundary layer 
and local chemical equi- 
librium at the vehicle 
surface. (See, for 

instance, ref. 1.) This 500 1000 1500 2000 2500 3000 3500 4000 °k 

assumption is equivalent iboo 2700 3 eoo 4500 5400 63 00 *r 

, . , , . Temperature , T 

to assuming a thin reac- Figure 3.- Range of validity for approximate solution, 

tion zone near the char 

surface. When reactions occur throughout the boundary layer, the present approach is 
not applicable; however, to the author's knowledge, no solutions for the transpiration 
effect of pyrolysis gases injected into a reacting boundary layer are available. Cer- 
tainly, there exists the possibility that transpiration cooling effects significantly different 
from those predicted by equations (120) to (127) will be obtained. As the following dis- 
cussion develops, however, it will be seen that transpiration factors which differed from 
the presently used values by factors of two or even three would not significantly change 
the qualitative results of this analysis. Hence the present assumption is reasonably 
accurate for the case of frozen boundary layers with chemical equilibrium at the char 
surface and could in no case introduce errors which would significantly alter the con- 
ceptual picture of diffusion-controlled char oxidation presented herein. 

The pyrolysis gases assumed in the present treatment represent the major species 
which have most often been reported by chemists who have attempted experimentally to 
determine pyrolysis gas compositions. (See, for instance, refs. 13 and 14.) The pyroly- 
sis gases considered in the present report are those which reach the char surface. 

Hence, while the initial decomposition of the plastic may produce large molecular weight 
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hydrocarbon species, it is assumed that cracking reactions in the char reduce these 
large fragments to the species listed in assumption (f). 

Pyrolysis gas compositions, calculated for the equilibrium decomposition (i.e., all 
cracking reactions in chemical equilibrium) of phenolic nylon are presented in refer- 
ence 15 and contain several species not included in assumption (f). The inclusion of all 
the species presented in reference 15 in the present analysis might possibly result in 
small changes in the region of validity shown in figure 3. The species included in 
assumption (f) do, however, account for approximately 90 percent (on a mass basis) of the 
gas mixture presented in reference 15 at temperatures near the middle of the range of 
validity presented in figure 3. At the low and high temperature extremes of the range of 
validity, the respective percentages of this gas mixture accounted for are approximately 
80 and 70. Hence, the presently assumed pyrolysis gases are considered to be a reason- 
able approximation for the case where the cracking reactions approach equilibrium. 

If, on the other hand, the cracking reactions were kinetically slow so that large molecular 
weight hydrocarbon species reached the char surface, the pyrolysis gases in assumption 
(f) would be incorrect and the presence of heavy gases would cause the boundary-layer 
gas mixture to deviate still further from an effectively binary mixture. The fact that the 
major pyrolysis gas species measured experimentally (refs. 13 and 14) and calculated 
(ref. 15) are the same strongly suggests that char layer cracking reactions approach 
equilibrium. 


Physical Interpretation of Approximate Solution 

A physical interpretation of the approximate solution may be obtained from an 
examination of equations (128) and (129). The numerator of equation (128) is the well- 
known expression for the diffusion-controlled oxygen flux and is the maximum rate at 
which oxygen can be transported across the boundary layer to the reacting char surface. 
Since equation (128) is an expression for the steady -state ablation rate, fW* must rep- 
resent the mass of oxygen used in ablating a unit mass of virgin plastic. If every carbon 
atom initially present in the virgin plastic is assumed to react with the oxygen in the vir- 
gin plastic to form CO, then the term — ^ represents the number 

crM^, + tM h + + (,c>Mq 

of moles of carbon per unit mass of virgin plastic available to react with oxygen diffusing 
through the boundary layer. Hence fW* (eq. (129)) represents the mass of oxygen from 
external sources required to react all the carbon in a unit mass of virgin plastic to pro- 
duce CO. This, then, is the physical interpretation of the approximate solution: all the 

carbon present in the virgin plastic is reacted to form CO and the reaction rate is deter- 
mined by the rate at which oxygen can be supplied by pyrolysis of the virgin plastic plus 
diffusion through the boundary layer. It is in this respect that the present analysis differs 
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most significantly from previous analyses of char oxidation and carbon combustion. In 
these previous analyses, only the carbon present in the char layer is considered reactable. 

Effect of Virgin Plastic Composition on Boundary -Layer Composition 
at the Char Surface, Char Oxidation Rate, Ablation Rate, 
and Effective Heat of Ablation 

The effect of virgin plastic composition on pyrolysis rate is shown by equa- 
tions (128) and (129). However, the transfer coefficient h which appears in equa- 
tion (128) is dependent on the virgin plastic hydrogen content as shown in equation (135). 
When equation (135) is substituted into equation (128), the following relation is obtained: 


rn 


m vp = 


*Q,e h ° 


K P,e h 


f fW* + Kq e rj fW*+ 0.7(1 + 1301^ Kq f 


(139) 


Equation (139) indicates that m^p decreases with increasing fW* and increasing 
transpiration factor rj. Since fW* increases with increasing virgin plastic carbon 
content and decreasing oxygen, nitrogen, and hydrogen content and ?j increases with 
increasing virgin plastic hydrogen content, the question arises as to whether the lowest 
pyrolysis rate would be obtained from a material with a high carbon content or a high 
hydrogen content. In an attempt to answer this question, equation (139) was used to make 
calculations of m V p for a material which contained only carbon and hydrogen. (Eq. (129) 
shows that the presence 
of oxygen or nitrogen in 
the virgin plastic always 
produces higher pyroly- 
sis rates.) The results 
of these calculations are 
presented in figure 4 as 
plots of m V p/h° 
against for 

% e = 0.050 and 0.232 
and l - 1/3. For a 
specified flow environ- 
ment ^h° and Kq g j, 
the minimum value of 
riiyp is obtained from 
a material which is 



Corbon moss fraction of virgin plostic .m^ 

Figure 4.- Variation of dimensionless steady-state ablation rate with carbon mass fraction for 
plastic which contains only carbon and hydrogen. 
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composed entirely of carbon. Numerous investigators have impregnated porous carbon 
materials with hydrocarbon resins in the hope that the transpiration cooling effect pro- 
duced by the vaporization of the resin would reduce the mass loss rate of the material. 

The results of this study show that, within the framework of the assumption presented 
earlier, the effect of impregnating a porous carbon with hydrocarbon resins would be to 
increase rather than reduce the mass loss rate. 

In order to illustrate the effects of virgin plastic composition, calculations of char 
oxidation rate, ablation rate, and char surface boundary -layer composition were carried 
out for four typical hydrocarbon resins and graphite. The empirical formulas and 
associated parameters used in these calculations are presented in table I. The indices 
in the empirical formulas have been adjusted to give the numbers of carbon and hydrogen 
atoms per atom of oxygen and the materials are listed in order of decreasing carbon con- 
tent. As pointed out previously the parameter fW* represents the mass of oxygen 
required to ablate a unit mass of virgin plastic. (For graphite, no pyrolysis is involved 
and hence fW* represents the mass of oxygen required to oxidize a unit mass of 
graphite.) The values presented in table I show that fW* increases markedly as the 
relative carbon content of the material increases. The symbol mjj represents the 
mass of hydrogen per unit mass of virgin plastic and is related to the transpiration 
cooling factor r\ by equation (125). The significant point to be made with the values pre- 
sented in table I is that while m H differs significantly for the five materials considered, 

TABLE I 

EMPIRICAL FORMULAS AND ASSOCIATED PARAMETERS FOR GRAPHITE 
AND FOUR TYPICAL HYDROCARBON PLASTICS 


Material 

Empirical 

formula 

fW* 

m H 

(a) 

Graphite 

C 

1.333 

0 

0.700 

Epoxy 

resin 

C 6.5 H 7.34 °1 

.869 

.0725 

.874 

Phenolic 

resin 

C 5.16 H 5 °1 

.802 

.0602 

.849 

Polyester 

resin 

C 2.5 H 2 °1 

.500 

.0416 

.810 

Cellulose 

C 1 2 H, 0, 

.099 

.0617 

.851 


a In the calculation of rj (see eq. (125)), the exponent l was assumed 
to be 1/3. 
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the maximum variation in q is only 25 percent. This again indicates that large 
increases in the transpiration cooling effect will not be obtained by increasing the hydro- 
gen content of a plastic. 


Graphite C(s) 


By using the empirical formulas presented in table I and equations (130), (131), 
and (132), the boundary -layer compositions at the char surface during the steady abla- 
tion of cellulose, polyester, phenolic, and epoxy resins, and a completely carbonaceous 
material were calculated and are presented in figure 5 as variations of Ng, Hg, or CO 
mass fraction with free -stream oxygen 
mass fraction. The char surface 
boundary -layer compositions are com- 
parable for all the materials except 
cellulose. Cellulose, with its high oxygen 
and hydrogen contents, produces much 
larger quantities of CO and Hg than the 
other materials do. The curves pre- 
sented in figure 5 show that (with the 
exception of cellulose) the boundary - 
layer composition for the plastics is 
similar to that produced by graphite 
(i.e., « K co and j. One of 

the principal assumptions involved in the 
present analysis is that the unity Lewis 
number heat-mass transfer analogy is 
valid. It is generally agreed that the 
unity Lewis number assumption is rea- 
sonably valid for carbon oxidation. 

(See ref. 9.) The comparable boundary- 
layer composition shown in figure 5 for 
epoxy and phenolic resins and graphite 
suggests that the unity Lewis number 
assumption may be a reasonable approx- 
imation for the present case. 



oxygen mass fraction. 


The effect of virgin plastic composition on steady-state ablation rate is illustrated 
in figure 6, where ratios of ablation rate to Kq e h° are plotted against oxygen mass 
fraction, Kq . The data are presented in this manner since the product Kq e h° rep- 
resents the maximum rate at which oxygen can be transported by diffusion to the char 
surface and hence characterizes the flow environment. The ratio e h ° indi- 

cates the ablator's vulnerability to oxidation. The curve for graphite is included since it 
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represents (as demonstrated earlier) the min- 
imum ablation rate attainable by a hydrocarbon 
material. The curves presented in figure 6 
show that plastics exhibit a range of carbon 
content which causes significant differences in 
char oxidation rates and ablation rates. The 
ablation rates calculated for phenolic and 
epoxy resins (the two plastics most commonly 
used in ablative heat shields) are about 1.4 to 
1.6 times as large as the theoretical minimum 
rates (for graphite), and the values for poly- 
ester and cellulose exceed those for graphite 
by factors of from 2.7 to 3 and 5 to 12, 
respectively. From time to time it has been 
suggested, in the literature, that wood (which 
is primarily cellulose) might be an acceptable 
ablator. The calculations presented in fig- 
ure 6, however, indicate that cellulose would 
not be an acceptable material for flight 
regimes where steady ablation and diffusion- 
controlled char oxidation are obtained. Of 
the three resins, phenolic, epoxy, and polyester, the present results indicate that phenolic 
and epoxy would be the least vulnerable to oxidation and indeed these are the two resins 
most commonly used in ablative heat shields. Even though phenolic s and epoxy s are 
relatively efficient ablators, the curves presented in figure 6 predict that completely 
carbonaceous materials (such as prepyrolyzed resins) would have mass -loss rates that 
were from 60 to 70 percent of the rates obtainable from phenolics and epoxys. Most 
carbonaceous materials have relatively high thermal conductivities and densities and 
hence are poor insulators. The proper interpretation of the curves presented in fig- 
ure 6 is, then, that if completely carbonaceous materials with insulating capabilities 
comparable to those of plastics could be developed, then these carbonaceous materials 
would be superior ablators by virtue of these lower mass loss rates. 

Estimates of the weight of material ablated by a given aerodynamic heat load can be 
obtained by defining an "effective heat of ablation." For present purposes a cold-wall 
effective heat of ablation is used and is defined as follows: 
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Figure 6.- Effect of virgin plastic composition on the ratio of 
ablation rate (or char oxidation rate) to diffusion oxygen 
flux. 
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m vp m vp 
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Substitution of equation (139) into equation (140) yields 


9e “ = (€r , ’) Hs <“*i> 

For diffusion-controlled oxidation, the response of the ablator is determined by the flux 
of oxygen to the char surface and is completely independent of the aerodynamic heat flux. 
As a result, equation (141) contains no thermal properties such as specific heats, thermal 
conductivities, or heats of pyrolysis. It is only in mass transport controlled cases, such 
as the present one, that an effective heat of ablation such as that defined in equation (140) 
has meaning. Where the response of the ablator is determined by aerodynamic heat flux, 
a thermochemical heat of ablation such as 



m vp 


(142) 


( where q net is the heat flux to a nonablating calorimeter that has the same tempera- 
ture, emissivity, and catalytic efficiency as the ablator ) must be used. 


Equation (141) was used to calculate values of Q e ^ for each of the materials 
appearing in table I for Kq = 0.232, which corresponds to air. The variation of these 
effective heats of ablation with stagnation enthalpy is presented in figure 7, which shows 


— Completely carbonaceous (graphite) 

Btu Epoxy 

lb Phenolic 



Stream stagnation enthalpy 


Figure 7.- Effect of virgin plastic composition on cold-wall effective heat of ablation. 
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that the cold-wall effective heats of ablation of epoxy, phenolic, polyester, and cellulose 
are, respectively, 72, 67, 45, and 19 percent of that for a completely carbonaceous mate- 
rial. Since, for a given cold-wall heat load, the weight of material ablated is inversely 
proportional to H e ^, figure 7 may be interpreted as showing that the respective weights 
of material ablated from heat shields composed of epoxy, phenolic, polyester, and cellu- 
lose would be 1.39, 1.49, 2.22, and 5.25 times the weight of material ablated from a com- 
pletely carbonaceous shield. 

Effect of Pyrolysis Gas Composition on the 
Heat Release Due to Chemical Reaction 

In order to illustrate the effects of varying pyrolysis gas composition, calculations 
were made for two assumed pyrolysis gas compositions. These compositions, which are 
presented in table II, were based on data presented in references 13 and 14. The results 

TABLE II 

ASSUMED PYROLYSIS GAS COMPOSITIONS 


Item 

Composition 1 
(phenolic resin) 

Composition 2 
(phenolic nylon) 

X »2 

0.5199 

0.2061 

*co 

0.1954 

0.2875 

x Na ...... 

0.0013 

0 

x co 2 

0 

0.0222 

^CH 4 

0.1427 

0.1980 

x c 2 H 2 

0 

0.1340 

Xc 2 H 4 

0 

0.0006 

x nh 3 

0.1407 

0.1515 

x h 2 o 

0 

0 

fW* 

0.724 

0.601 

m H 

0.0974 

0.0795 

7 

0.921 

0.886 


C 5.9 H 10.37 N 0.73 °1 

C 3.74 H 5.80 n 0.46 °1 

f 

0.462 

0.23 
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of these calculations are compared with those for graphite and are presented in fig- 
ures 8, 9, and 10. 

From table II it is seen that the two compositions differ significantly with respect 
to Hg, CO, and CgHg content. A comparison of table II with table I shows that the values 
of fW* and the virgin plastic carbon and oxygen contents derived from the pyrolysis gas 
compositions are lower than those derived from the molecular structure of phenolic 
resin. Because of the extreme difficulty of experimentally determining pyrolysis gas 
composition, and because, in the compositions presented in table II trace species and 
condensibles included in references 13 and 14 were neglected, the compositions presented 
in table II must be considered as assumed rather than experimentally measured. Con- 


sidered in this light, the agreement in fW 

Figure 8 presents the boundary - 
layer composition at char surface as a 
function of Kq for two assumed com- 
positions and graphite. Comparing fig- 
ure 8 with figure 5 shows that the char 
surface boundary-layer composition pro- 
duced by the assumed pyrolysis gases 
closely approximates that derived from 
the molecular structure of phenolic 
resin. 

In figure 9 the variation of the 
maximum heat of reaction (per pound of 
char) with char surface -temperature is 
presented for graphite and the two pyrol- 
ysis gas compositions. The maximum 
heat of reaction AH^ for composition 1 
is only slightly lower than that for graph- 
ite, whereas AHj^ for composition 2 is 
approximately twice that for graphite. A 
large part of this difference is due to the 
difference in f (mass of char produced 
when a unit mass of virgin material is 
pyrolyzed) between compositions 1 and 2. 
(See eq. (116).) Figure 9 shows that AH^ 
(per pound of char) can be a strong function 
of virgin plastic composition and pyrolysis 
gas composition. 


is perhaps not too bad. 



Figure Variation of ooundary- layer composition at char surface 
with free-stream oxygen mass fraction as calculated from 
assumed pyrolysis gas compositions. 
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Figure 9.- Variation of maximum heat of reaction with char surface temperature as calculated from assumed 

pyrolysis gas compositions. 


As shown by equations (117), (118), (119), and (134) an effective enthalpy potential 
due to chemical reaction can be defined as a function of virgin plastic composition, pyrol- 
ysis gas composition, and free -stream oxygen mass fraction. This effective enthalpy 
potential ( per unit Kq \ is plotted against char surface temperature in figure 10. Note 



1—1 1 1 I I _L i J I J 1 — J I 1 I J... i 1 I J 


0 1000 2000 3000 4000°K 

I I 1 I J _ 1 I 

1000 2000 3000 4000 5000 6000 7000 °H 


Char surface temperature , T w 

Figure 10.- Variation of effective enthalpy potential due to oxidation with char surface temperature as 
calculated from assumed pyrolysis gas compositions. 
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that AH c j i for composition 2 is larger than that for graphite whereas AH^ for com- 
position 1 is less than that for graphite. The parameters responsible for these differ- 
ences can be identified by examining equation (134). Rearrangement of equation (134) 
yields 






1,00 ' x co + 2X C0 2 * x h 2 o 


where the first term on the right-hand side is the effective enthalpy potential for the com- 
bustion of graphite. Note that the second term on the right-hand side will be positive or 
negative ^and hence AH ch will be greater than or less than ^AH ch ^^j depending on 

l VV 

whether AH, nr . is greater than or less than — i-rr — . Also note that 

Y j- 2Y a. Y__ _ 


x co + 2 x co 2 + x h 2 o 


l V H M 

for a given value of AH, = the magnitude of AH , depends 

Y « OY . X LI1 


co T co 2 T "h 2 o 


on the size of the ratio ^ . Since AH f ^ is a negative quantity, AH ch will be 

larger than ^AH ch ^ ^ when ^ XAHj i is positive. Of all the pyrolysis gas con- 
stituents considered in this paper, only C 2 H 2 and CgH^ have positive values of AH. 

Hence AH ch becomes larger as the amounts of C 2 H 2 and C 2 H 4 in the pyrolysis gases 
increase. Note from figure 10 and table H that AH ch for pyrolysis gas composition 2 
(which contains C 2 H 2 ) is larger than AH ch for pyrolysis gas composition 1 (which 
does not contain CgHgJ. Vojvodich and Pope reported in reference 17 heating due to 
chemical reactions much larger than could be accounted for by carbon combustion. From 
equation (143) it is seen that very large values of AH^ could be produced by a plastic 
which had a small value of (a - q>) and produced pyrolysis gases containing large amounts 
of C 2 H 2 and C 2 H 4 . 
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SUMMARY OF RESULTS 


\ 


An analytical study of diffusion-controlled char oxidation and its effect on the abla- 
tion of hydrocarbon plastics has been carried out. Although the solution obtained is 
based on a number of simplifying assumptions, it is believed to present an accurate con- 
ceptual picture of the diffusion -controlled char -oxidation process. The steady -state 
expressions resulting from this analysis can be used directly to provide a clear physical 
picture of char oxidation. This analysis has also produced general expressions which 
can be used in conjunction with a transient heat conduction computer program to yield 
transient ablation and char oxidation rates, but such a procedure does not lend itself to 
demonstrating the general characteristics of diffusion -control led char oxidation. Accord- 
ingly, the discussion and results presented in the present report were restricted to the 
steady-state case. The results of this study are summarized as follows: 

1. The present analysis differed from previous analyses of diffusion -controlled 
char oxidation in one important respect: previous investigators have assumed that only 
the carbon in the char layer would react and that the rate of reaction would be determined 
by the diffusion rate of oxygen through the boundary layer but the present analysis has 
shown that the carbon present in both the char and the pyrolysis gases (i.e., all the 
carbon originally present in the virgin plastic) reacts and that the reaction rate is 
determined by the rate at which oxygen can be supplied by diffusion through the boundary 
layer plus pyrolysis of the virgin plastic. Because of this fact, the char oxidation rate, 
the steady-state ablation rate, the boundary -layer composition at the char surface, and 
the heat released by chemical reactions were found to be dependent on the composition 

of the virgin plastic and the pyrolysis gases. 

2. For a given flow environment (i.e., oxygen mass fraction and heat -transfer coef- 
ficient) the char oxidation rate and the pyrolysis rate were found to increase as the 
hydrogen, nitrogen, and oxygen contents of the virgin plastic increased and to decrease 
as the carbon content of the virgin plastic increased. Even when transpiration cooling 
effects were accounted for, the lowest steady -state pyrolysis rate was obtained from a 
material composed entirely of carbon. 

3. The cold-wall effective heat of ablation of a plastic material was found to be 
determined by the free-stream oxygen mass fraction, the stream enthalpy potential, and 
relative number of carbon, hydrogen, nitrogen, and oxygen atoms in the virgin plastic. 
Example calculations carried out for a completely carbonaceous material, epoxy, phenolic 
and polyester resins, and cellulose showed that the effective heats of ablation for the 
polymers would be, respectively, 72, 67, 45, and 19 percent of the effective heat of abla- 
tion calculated for the carbonaceous material. 
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4. The maximum heat of reaction per pound of char was a strong function of 
pyrolysis gas composition and the fraction of virgin material which was converted to 
gas, and the heat released by both gas phase and heterogeneous reactions can be 
accounted for by an "effective enthalpy potential" which is a function of the stream 
oxygen mass fraction, the relative number of carbon, oxygen, hydrogen, and nitrogen 
atoms in the virgin plastic and the composition of the pyrolysis gases. The pyrolysis 
gases considered in this investigation were CO, COg, H 2 , N 2 , NHg, HgO, CH 4 , C 2 H 2 , 
and CgH^. Both the maximum heat of reaction and the effective enthalpy potential were 
found to increase with increasing C 2 H 2 and C 2 H 4 content. Particularly high heat 
releases were predicted for polymers which contained nearly equal numbers of carbon 
and oxygen atoms and produced pyrolysis gases containing large quantities of C 2 H 2 and 
and C 2 H 4 . 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., February 18, 1966. 
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